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We demonstrate the advantages of randomization in coherent quantum dynamical control. For 
systems which are either time-varying or require decoupling cycles involving a large number of 
operations, we find that simple randomized protocols offer superior convergence and stability as 
compared to deterministic counterparts. In addition, we show how randomization may allow to 
outperform purely deterministic schemes at long times, including combinatorial and concatenated 
methods. General criteria for optimally interpolating between deterministic and stochastic design are 
proposed and illustrated in explicit decoupling scenarios relevant to quantum information storage. 
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Counterintuitive roles of randomness in physical pro- 
cesses have long been acknowledged. A paradigmatic ex- 
ample is stochastic resonance 0], whereby the transmis- 
sion of a weak signal is maximized at an optimum noise 
level. Within quantum information science, a number 
of applications benefiting from randomness in quantum 
states and operations have been identified recently. Sug- 
gestive evidence is offered bv dissipation-assisted quan- 
tum computation schemes ^ and by the idea that chaos 
may stabilize quantum algorithms 3J. In the quantum 
communication setting, the use of random operations 
has been shown to decrease the communication cost of 
achieving remote state preparation and of constructing 
efficient quantum data- hiding schemes 0|. Finally, ran- 
dom unitary operators may pave the way to efficient es- 
timation methods for open quantum systems f^. 

In the context of quantum control theory, the fragility 
of quantum coherence against uncontrolled noise and 
randomness hindered the exploration of stochastic con- 
trol methods till recently |E El B Q- General 
control-theoretic arguments supporting the usefulness of 
controlled randomness in dynamical decoupling (DD) 
schemes for switching off unwanted evolution in interact- 
ing quantum systems were provided in ^] , and have been 
validated in specific examples |^ (see also l9]). These 
studies leave, however, several key questions largely 
unanswered. In particular: For given control resources, 
under what circumstances (if any) can randomized DD 
match or outperform the best deterministic schemes 
known to date? Can randomized design be exploited in 
general to further push the efficiency and viability of DD 
schemes? 

In this Letter, we present conclusive evidence of the 
benefits of randomization in light of the above questions. 
To isolate the factors responsible for different DD per- 
formance, we focus on the simplest control benchmark: 
complete refocusing of Hamiltonian evolution by means 
of a restricted but otherwise perfect set of control oper- 
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ations. In the process, we verify conjectures made in [3], 
stretch the analysis into unexplored domains, and offer 
criteria for constructing new highly efficient DD proto- 
cols. For time-independent systems, we demonstrate the 
superior averaging of randomized protocols in control sce- 
narios involving sufficiently long control cycles and/or 
evolution times. For time- varying systems, randomized 
DD provides a robust strategy in the presence of system 
uncertainties. The combination of three key ingredients 
- concatenated control, symmetry, and randomization - 
emerges as a general principle for design optimization. 

DD setting.- DD methods modify a target dynamics 
Ho{t) through the addition of a control field Hc{t). We 
assume that ||iJo(i)||2 = max \eig{Ho{t))\ < k, for k > 
and all t. Let Uc{t) — Texp[— i Hc{u)du] {h = I) 
denote the control propagator, with T indicating time 
ordering. In a logical frame which explicitly removes 
Hc{t), the controlled evolution is described by the prop- 
agator U(t) = Ul{t)U{t) = T exp[-i f^H{u)du], where 
U{t) is the total propagator in the physical frame and 
H{t) = Ul{t)Ho{t)Uc{t) is the logical Hamiltonian lOj. 
If Ho{t) is time-independent and Uc{t) periodic with a 
cycle time T^, then U{Tn) U{Tn), r„ = nT^, n e N, 
and both frames coincide stroboscopically. The DD ob- 
jective here is to make U{T), for final time T > 0, as 
close as possible to 1 according to a metric of choice. 

Let the control resources be specified by a discrete set 
of unitary operators corresponding to a (projective) rep- 
resentation of a group G = {gj}, J = 0, . . . , |^?| — I, in the 
Hilbert space H of the system, dim(7Y) ~ d < 00. Both 
the basic deterministic and random DD protocols may be 
understood as effecting an appropriate symmetrization 
of Holt) according to Q. Assume a time-independent 
system first. In standard deterministic DD, average 
Hamiltonian theory (AHT) allows the logical evolution 
to be represented as i7(T„) = [exp(— i/fTc)]", where 
H =X]fe°=o ^iid each term H^^' is computed from the 
Magnus expansion [lol | . Under the convergence condition 
kTc < 1, the leading zero-th order contribution H^-"^ may 
be mapped into a group average Hg = \G\~^ dj^agj 
by sequentially steering Uc{t) through all the group ele- 
ments {gj} - which corresponds to a train of "bang-bang" 
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instantaneous pulses Pk = gkgl-ij ^ = 1; ■ ■ ■ j sepa- 
rated by A< > and with ^ \g\At A cyclic pro- 
tocol based on a fixed, predetermined control path within 
Q will be referred to as periodic DD, henceforth (PDD), 
provided first-order averaging is achieved, = 0. The 
simplest randomized version of such a DD procedure is 
obtained by sampling Uc{t) over Q uniformly at random 
according to the Haar measure Q - leading to what we 
call naive random DD (NRD). Let E{-} denote expectation 
over all possible control realizations; Cy-symmetrization 
still emerges on average at each t via a quantum opera- 
tion E{U^{t)Ho{t)Uc{t)} ~ Hg, and convergence may be 
rigorously established in the limit TAtu^ <C 1 [3 • 

Lower bounds for the expected logical-frame fi- 
dehty of an arbitrary pure state \t{j), E{F|^)(T)} = 
E{|(V'|t/(T)|^}|2}, are given by E{F|^)(T)} > 1 - 
0{T^T^K*) for PDD, and E{F|^)(r)} > 1 - 0{TAtK^) 
for NRD, respectively 0. Within their regime of validity, 
these bounds indicate the potential for NRD to outperform 
PDD when \Q\'^(TAtK'^) > 1. In order to quantitatively 
compare DD schemes, a control metric which is both ef- 
ficiently computable and state-independent is desirable. 
Here, we consider average (unlike worst-case as in 0) 
performance, and remove the dependence upon lip) by 
invoking the (expected) entanglement fidelity 12], evalu- 
ated as E{Fe(r)} = E{|tr(l7(T))/dp}. FeiT) is linearly 
related to the average of (T) over all Perfect DD 
corresponds in this metric to E{Fe{T)} ^ 1. In practice, 
Monte Carlo simulations of the controlled dynamics are 
used to estimate E{Fe{T)} through a statistical average 
((Fe(T))) over control realizations. 

Time-independent case: Convergence improvement- 
By way of illustration, consider N spin-1/2 particles 
(qubits), described by the following model Hamiltonian: 

N N 
2—1 i<j a—x,y,z 

where cr^^^V'^^ = X,Y, Z are Pauli operators, and Ui, Jij\ 
denote the frequency of the «th qubit, and the cou- 
pling strength of the ijth pair in the a direction, re- 
spectively. Heisenberg interactions exponentially decay- 
ing with the distance rij are typical of quantum dot ar- 
rays ; whereas cubic decays describe dipolarly coupled 
spins in systems ranging from nuclear magnetic resonance 
(NMR) crystals and liquid- crystals 0,13 to electrons 
on Helium Ts^. We assume here that uii ^ uj, and work 
in a logical-rotating frame whereby the effective Hamilto- 
nian Hnit) = U}it)ulit)[H{t) - wY.,Z,/2]UR{t)Uc{t), 
UB.{t) ^ exp[~iu;tY.^ ZJ2] ^E^. 

PDD protocols capable of refocusing H for arbitrary pa- 
rameter values (and in fact, arbitrary interactions of the 
form cr|°^ (g) uj''^ ) may be built by recursively nesting DD 
sequences based on the group Qi = (l,; , X j, Zj, Yj} for 
each added qubit, i = 2,...,N U^M,U^. Although 
this scheme is not efficient, as the number of pulses per 
cycle grows as 4^~^, it allows the effect of large \Q\ to 




FIG. 1: (color online) PDD (dashed lines) vs. NRD (solid line) 
based on a nested pulse sequence for H from Eq. Q in ID 
with A*' = 6, J^j"' = J\i ~ in the logical-rotating frame. 
Main panel: average fidelity at T„ = n\Q\At. Free evolution: 
oscillating solid line. Inset: average fidelity within a cycle, 
tn ~ nAt, At — 10~'^J~^. Dotted line: average over different 
path choices (RPD, see text). Average over lO'^ realizations. 

be studied in moderately small systems. Following [l8l |. 
the PDD sequence we implement corresponds to a path 
over the Pauli group = ®iQi which avoids simultane- 
ous pulses. In NRD, Uc{t) is picked uniformly at random 
over , which may involve collective rotations on up to 
R ^ N — 1. qubits j^^l . Numerical results are shown in 
Fig. n As seen in the main panel, the fidelity for NRD de- 
cays substantially slower than for PDD. Irrespective of the 
validity of the short-time condition underlying the lower 
bounds of 7] , these results confirm the faster convergence 
of stochastic DD when |5| is large. While At is kept fixed 
for NRD, PDD protocols with decreasing At are considered. 
NRD eventually surpasses all PDD curves, showing that, for 
sufficiently long time, the constraints on At for random 
DD may be relaxed - which can prove advantageous in 
practice. Also note that in situations where control con- 
straints make it unfeasible to complete a cycle (inset), 
NRD outperforms the selected PDD sequence at all intra- 
cycle times t„ = nAt < Tc - a finding which indicates 
the existence of better control paths in this time domain. 

Time-independent case: Long-time improvement- 
While different DD protocols are equivalent in the short- 
time limit under ideal averaging, performance for finite 
At and long time depends critically on the way in which 
residual errors accumulate. We address this point by ex- 
tending our comparison to higher-level DD schemes. 

In the deterministic domain, two strategies exist for 
improving over PDD. One, motivated by the Carr-Purcell 
sequence of NMR, consists of a time-symmetrization of 
the control p ath, such that all odd order terms in H are 
also canceled ^JEl- This leads to symmetric DD (SDD). 
The other, inspired to NMR iterative techniques 10], im- 
plements concatenated DD (CDD) |^. CDD relies on a 
temporal recursive structure, the {i H- l)-level pulse se- 
quence being determined by Cf+i = CePiCiP2 ■ ■ ■ CiPn, 
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FIG. 2: (color online) Deterministic vs. randomized DD 
based on Qs for Hamiltonian H as in Fig. except A*' = 
\Q\ = 8. Data are averaged over 10^ control realizations, 
At = 0.05J~^. Note that at times r„ — 8nAt, as considered, 
the inner sequences characterizing SDD, CDD, and SRPD are not 
necessarily completed. For CDD, ^ = 5 is achieved at 102. 4J^^ . 
Dashed and dot-dashed lines correspond to different inner PDD 
sequences for EMD. Free evolution: (black) solid line. 

where Pk is the fcth pulse, Co is the inter-pulse interval, 
and Ci denotes the generating PDD sequence. 

Deterministic DD is plagued by coherent error build-up 
at long time - as opposed to probabilistic error accumula- 
tion in NRD. In order to construct and classify stochastic 
DD schemes able to ensure good performance at both 
short and long times, we describe each protocol in terms 
of an mner and an outer code. The former determines the 
control path in Q, thereby the pulse sequence in the inter- 
val [n, {n + \)\\Q\/S.t, subject to the condition that a net 
effective Hamiltonian with norm 0(At) is obtained. The 
outer code specifies additional pulses to be applied at T„ 
according to an outer group Q' , which need not coincide 
with Q. Combining a fixed PDD sequence with border- 
ing pulses drawn at random from Q' leads to embedded 
DD (EMD), as in with g' = . While choosing g' 
irreducible typically ensures superior performance, EMD 
schemes based on reducible g' — g may still prove useful 
under restricted control capabilities. However, a poten- 
tial weakness of EMD is non-uniform performance across 
the set of \g\\ inner paths. Path optimization may be 
shortcut if, at every r„, a random path choice is effected 
- resulting in what we term random path DD (RPD) j^.l22j|. 
Randomization of the inner code ensures that DD perfor- 
mance are robust against path variations. To further im- 
prove, control paths which are symmetrized in the same 
manner as in SDD may be employed, leadiiig to a scheme 
we call symmetric random path DD (SRPD) [23|. As a final 
step, interpolated DD protocols, which smoothly switch 
from high-level deterministic DD (e.g. CDD) to high-level 
randomized DD (SRPD), allow to optimize control perfor- 
mance over the entire time axis. 

Qualitative insight into the above-mentioned DD pro- 
tocols may be gained through analytical bounds on the 
expected fidelity decay. For deterministic schemes, we 



may write E{Fe{T)} > 1 ~ 0{T'^(\g\At)^+" k'^+°'), where 
a = 0, 2 for PDD Jj and SDD JJi], respectively. Error 
bounds for CDD depend on both the system and the con- 
catenation level, but already at £ = 2 CDD is more efficient 
than SDD in reducing H'^''\ k >2. Notice, however, that 
unlike the single-qubit DD analyzed in 21], no guarantee 
of superpolynomial convergence exists when CDD is based 
on a reducible group. For randomized DD, the basic error 
bound is contained in for NRD, further adapted in Q 
to EMD. While full derivations will be presented in we 
find E{Fe(T)} > 1 - 0(T(|g|At)3+"«;4+a)^ ^Yiere a = 
for EMD and RPD, while a = 2 for SRPD. As applicability of 
the above bounds is strictly confined to the regime where 
T(|CJ|A<)k^ ^ 1, we proceed to a numerical comparison. 

For bilinear couplings as in Eq. efficient combi- 
natorial PDD schemes with quadratic complexity are well 
known 18]. First-order DD sequences for up to iV = 4'" 
qubits, m > 1, involve N simultaneous pulses. For 
N = 8, the basic PDD scheme is derived from the group, 

Z2Y3X4ZQYJXS, Z2X3Y4Yr-,XQZj, Y2YiX^Z^X-;Z^, 

i^2^3^4^5-'^6^8i -'^2^3^4-^5^7^8i -'^2-'^3^5^6^7^8}- 

As C/8 is reducible, saturation in the performance of CDD is 
observed when £ > 2. Numerical results are summarized 
in Fig. [3 As expected, randomized DD surpasses deter- 
ministic DD at long times, the crossing being evident be- 
tween protocols that have comparable efficiency at short 
times. The dependence of EMD upon the underlying inner 
path is very pronounced for this system. Remarkably, 
CDD is outperformed by SRPD. Although CDD is known to 
minimize sensitivity to control faults not included here, 
preliminary results on rotation angle imperfections indi- 
cate that the above conclusions remain unchanged. 

Exploitation of randomized DD in actual devices will 
vary depending on whether probabilistic pulse generation 
capabilities are available or not. In the former case, ran- 
domization appears especially promising wherever a large 
number of experimental runs is required for observation, 
such as in quantum dot ensemble experiments 12^ ] . Even 
in the absence of dedicated pulsing capabilities, random- 
ized DD may still be useful (e.g., in solid-state NMR 
systems) once the procedure is de-randomized^ that is, 
randomization is used off-line to post-select optimal pulse 
realizations for given physical parameters. Such an op- 
tion is no longer viable for time- varying control systems. 

Time-varying systems and robust performance.- For 
simplicity, let us specialize Eq. ^ to a situation 
where only nearest-neighbor couplings are relevant, 
but a time-dependent anisotropy is effectively present, 
J^fj^i = J, Jf^i+i — JA{t). For arbitrary N, an 
efficient first-order PDD protocol may be constructed 
by alternating two collective rotations around per- 
pendicular axes - one acting on all odd qubits, the 
other on the even ones, e.g., if N is even, g — 
{1, Z1Z3 . . . Z^^i, Z1Y2Z3Y4 . . . Zjv-iIat, Y2Y4 . . . Y^v}. 

While under mild assumptions on the time depen- 
dence of A(t), the NRD fidefity bound is stiU vafid 0], 
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FIG. 3: (color online) Deterministic vs. randomized DD 
based on a —4 sequence for ID nearest-neighbor 

couplings with A'' = 8 in the logical-rotating frame. 
A{t) =X]fe=i sm{10TvRkJt), Rk uniformly random in [0.9, 1.1]. 
Average l^delity at T„ = n\Q\At, 10^ realizations. PDD: 
down-triangles; SDD: up-triangles, CDD: squares; NRD: solid- 
line; RPD: long-dashed line; SRPD: dot-dashed line. Free evolu- 
tion: (black) short-dashed line. Main panel: At — 0.05J~^, 
t/ At ~ 2. Inset: At = 0.025J~^ - notice the sensitivity of 
deterministic protocols (including CDD) to changes in At. 

no general result on deterministic performance is avail- 



able. Suppose that the system fluctuation entails sign 
changes over time, with a potentially unknown period 
T. Then if r/At ^ p, p £ N, adversarial situations may 
arise where a pre-established control action is inhibited or 
even reversed, resulting in surprisingly poor performance 
of cyclic DD - often worse than NRD or free evolution, see 
main panel of Fig.|3| Stochastic methods are intrinsically 
more protected, on average, against such interferences, 
resulting in more stable DD performance throughout. 

Conclusion - We have shown that randomized decou- 
pling can offer distinctive advantages over deterministic 
methods in terms of faster convergence, long-time im- 
provement, and robust performance. Together with con- 
catenation and symmetry, randomization provides a ver- 
satile toolbox for matching different control needs. Both 
the analysis of fault-tolerance properties and of the inter- 
play between probabilistic design and robustness deserve 
a closer scrutiny, also in view of known randomized al- 
gorithms for classical control systems j25| . It is our hope 
that the present results will prompt experimental verifi- 
cation in available control devices. 
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